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Abstract 

We study the Sivers function for valence u and d quarks in the proton within the in- 
stanton model for QCD vacuum, adopting the MIT bag model wave functions for quarks. 
Within approaches based on perturbative one-gluon final state interactions a non-zero 
value of the Sivers function is related to the presence of both S and P wave components 
in quark wave functions. We show that the instanton-induced chromomagnetic, non- 
perturbative interaction leads to very specific spin-spin correlations between the struck 
and spectator quarks, resulting in a non-trivial flavour dependence of the Sivers function. 
Comparison of the obtained Sivers functions with phenomenological parameterizations is 
discussed. 
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1 Introduction 



Understanding the spin structure of the nucleon is a long-standing, challenging problem 
for strong interaction theories which has received renewed interest in recent years. The 
investigation of mechanisms able to explain the unambiguous observation of large trans- 
verse single spin asymmetries (SSA) in inclusive meson production in hadronic collisions 
and in semi-inclusive deep inelastic scattering (SIDIS) may help in understanding the 
origin of spin effects in hadronic physics. Recently, significant SSA's in pion electro- 
production on both longitudinally [1], and transversely [2] polarized proton targets have 
been observed by the HERMES Collaboration at DESY. Analogous measurements and 
processes are actively under investigation by the COMPASS experiment at CERN [3]. 
Measurements of the SSA for pion production in polarized proton-proton collisions, both 
in the central rapidity region and at large Feynman xf, moderately large transverse mo- 
mentum and center of mass energy ^/s = 200 GeV, have been performed by the PHENIX 
[4], STAR [5] and BRAHMS [6] Collaborations at the Relativistic Heavy Ion Collider 
(RHIC) at Brookhaven. These experiments plan to refine their measurements, improving 
the statistics and enlarging the kinematical region explored. 

We should emphasize that in spite of the fact that large, "anomalous", SSA at high 
energy and large transverse momentum were discovered long time ago in both inclusive and 
exclusive hadron production in hadron-hadron collisions [7], their understanding within 
QCD [8, 9, 10], although significantly improved in recent years [11, 12], is far from being 
complete, see e.g. Ref. [13]. 

In perturbative QCD (pQCD) approaches with inclusion of intrinsic (or primordial) 
parton motion, SSA can be related to the so-called naively T-odd, transverse momentum 
dependent (TMD), Sivers distribution and Collins fragmentation functions [14, 15, 16]. 
These functions provide very important information on the behaviour of strong interac- 
tions at large distances and can be calculated only within nonperturbative QCD. There- 
fore, in the absence of first-principles QCD calculations for these functions, e.g., within 
lattice QCD, it is very important to calculate SSA within some suitable model accounting 
for nonperturbative QCD effects. In the last years it has been shown that the instanton 
model (IM) can well describe the main nonperturbative aspects of QCD [17, 18, 19, 20, 21]. 
The first attempts to apply the IM to the study of SSA were done in Ref.s [22, 23, 24, 25]. 
In particular, the instanton induced gluon-quark interaction was shown to be important 
for the appearance of significant SSA in quark-quark scattering [22, 23]. In Ref. [25] this 
mechanism has been applied to estimate the current jet single spin asymmetry in SIDIS, 
within a modified version of the final-state interaction approach suggested in Ref. [11]. 
The contribution to the SSA in SIDIS coming from the instanton-induced spin-flip photon- 
quark vertex is the subject of a recent paper [24]. 

In this paper we derive for the first time the Sivers distributions for u and d valence 
quarks inside a proton, within the instanton model for QCD vacuum and using the MIT 
bag model for quark wave functions. Part of the calculation is performed along the lines of 
Ref. [26], where however the perturbative one-gluon exchange contribution to the Sivers 
function was considered. Analogous calculations, using a quark-diquark model for the 
nucleon, were presented in Ref.s [27, 28]. 
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2 Sivers function within the MIT bag model formal- 
ism 

The Sivers function [8], fi T a (x, k\) or A N f a/p1 (x, k 2 ± ) (see Ref. [29] for details on differ- 
ent notations adopted in the literature), describes the correlation between the intrinsic 
transverse momentum of an unpolarized quark (with flavour a) and the parent proton 
transverse polarization. The general correlator for transverse momentum dependent quark 
distribution functions is given by [14] 

* a (x, k ± , s) = J e -*< (P , S \ r(ov\ovm a (o) |p, s)y =0 , (1) 

where k^ = xP^ + k^_ is the quark four-momentum, and P, S are the hadron momentum 
and spin vector respectively. The Wilson lines, also called gauge-links, V(£), are path- 
ordered exponentials required to preserve gauge invariance; as it was shown recently, they 
are also crucial in order to get non-vanishing T-odd TMD parton distributions, like the 
Sivers function [11, 12]. They are defined as follows: 

V(0 =Vexp(-ig s J™dz- A + (z-,^ ■ V exp (-ig s J™ dz ± A ± (z~ = oo, z ± )^) , (2) 

where g s is the strong coupling constant, = I] a t a A^, t a = A a /2 are the Gell-Mann 
colour matrices and A% are the gluon gauge fields. In the following, we will work in a 
covariant gauge, so that we can neglect the second term in Eq. (2), which vanishes in 
nonsingular gauges [30]. 

TMD distribution functions can be obtained by first opportunely Dirac-projecting and 
then taking the trace of the general correlator, Eq. (1). In the case of interest, for example, 
the Sivers function for a quark with flavour a can be extracted as follows: 

ft?(x, k ±) = o^ Tr 7 + k ± , S ±x ) - ® a (x, k ± , -S ±x )] , (3) 

where k± = \k±\, M P is the proton mass and, to be definite, we have considered a proton 
moving along the +z direction and transversely polarized along the x direction; in this 
way, we have also subtracted the fcj_-even contribution. Due to translational invariance, 
using Eq. (1) one can rewrite the last equation in the form 



(4) 

Hereafter, for shortness, we will indicate by ( \ x | • • • | \ x ) any matrix element of the 
form [ (S± x \ ■ ■ ■ \S± X ) - (S± x \ ■ ■ ■ \ - S ±x ) ]/2. 

In order to estimate the instanton contribution to the Sivers function, following Ref. [26] 
we will use (valence) quark wave functions from the MIT bag model. The Fourier trans- 
form of the lowest modes of these functions, to which we limit, is then defined as 

^(0 = E/S^)^< m , (5) 

where a l a m is the quark annihilation operator, a, i, are respectively the quark flavour and 
colour indexes and m is its spin component (helicity), m = ±1/2. Since we consider only 



2 



valence quarks in the proton, a similar term related to antiquark operators is omitted. 
Moreover, in the MIT bag model 



<f> m (k) = nV^r 3 



to(k) Xm 
a k ti(k) Xr 



(6) 



where k — \ k |, k — k/k, 

f 1 1 

U(k) = ^ dx x 2 ji(xkR )ji(xuj) , N 2 = r , ( , 



2(u — 1) sin ijj 



(7) 



i?o is the bag (confinement) radius, and uj = 2.043 for the lowest state of massless quarks 
inside the bag. In the proton rest frame Ep = Mp and the bag radius i?o = Auj/Ep is 
fixed by the bag stability condition, dE P /dR = 0. 

The transversely polarized proton ( e.g., in the ± x-direction) is described by the state 
vector 

(8) 



\S ±X ) = -j= (\S +Z ) ± , 



where the longitudinally polarized states read 

^U^iKlMalVs) - <tK)<I {n 2 )aj{m)\ | 0) , (9) 
~ a ui Ma J u \(n 2 )aJ(n 3 ) + a^n^a 3 } (712)0^(713)] I 0) , 



1 1 



y/E y/l8 L ™l,™2,n3 

and the rii are position indexes 



3 Instanton contribution 

The contribution to the quark Sivers function coming from perturbative, one-gluon ex- 
change, final state interactions (see Fig. la) has been considered in Ref. [26] in the frame- 
work of the MIT bag model. The interaction Lagrangian in this case is simply given 
by 

£L = -EE s.^yti^w^) ( 10 ) 

fS a,k,l 

where (5 and a, k, I are respectively flavour and colour indexes. It has been shown in 
Ref. [26] that the interference between S and P wave components of the bag quark wave 
functions leads to a non-zero value of the Sivers function. Clearly, these contributions 
are related to the dynamics of quark confinement at the scale of a typical hadronic size; 
as such, they can be responsible for the SSA at relatively small transverse momentum, 
k± « 1/R « 200 - 300 MeV/c. 

In this paper we consider the contribution to the Sivers function which stems from 
a completely different QCD dynamics (see Fig. lb), namely from the nonperturbative 
chromomagnetic quark-gluon interaction induced by instantons [22] 

££f = i — E E A ^' k (v) <v tli ^' l (v) <%r(v) , (ii) 

9s p a ,k,i m (3 
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Figure 1: Perturbative (a) and instanton induced (b) contributions to the Sivers function. 
The symbol I denotes the instanton. 

where / = ir 2 n c p^ is the so-called packing fraction of instantons in the QCD vacuum, 
n c is the instanton density, p c is the average instanton size, m* ~ m* d ss m* is some 
appropriate mass parameter for the quark propagator in the instanton vacuum [17] and 
a iiv = bVt,7^]/2. In mean field approximation for free quarks m* has the meaning of 
an effective quark mass (see, e.g., the recent discussion in Ref. [31]). For off-shell quarks 
with virtualities k\ and k 2 and a gluon with virtuality q the instanton induced quark-gluon 
vertex should be multiplied by the factor 

F(h, k 2 , q) = <S>(k lPc /2)<S>(k 2 p c /2)G(qp c ) , (12) 

where 

$(z) = -z-^[I (z)K (z) - IMK^z)] , 

G(z) = - 2 -1K 2 {z) (13) 

z z 

are the Fourier-transformed quark zero-mode and instanton fields respectively and I u (z), 
K u (z), are the well-known modified Bessel functions. The main feature of the instan- 
ton vertex, Eq. (11), in comparison with the perturbative one, Eq. (10), is the presence 
of a quark spin-flip contribution. Moreover, the strong coupling constant g s enters the 
denominator of the effective strength in Eq. (11); therefore there is not any a s power 
suppression of the spin-flip contributions in our nonperturbative approach to SSA. We 
should also mention that the average instanton size is much smaller than the confinement 
radius, p c ~ (l/3)-R co „/, so that the instanton contribution to the SSA should be sig- 
nificant at larger momentum transfer in comparison with the contribution coming from 
quark bound-state dynamics. 

Let us now estimate the contribution of the chromomagnetic interaction to the quark 
Sivers function coming from Fig. lb. Inserting Eq. (11) and the gauge links, Eq. (2), into 
Eq. (4), and taking the leading term we find 

,i / , x MpE P n 2 f [d£-d 2 £ ± _ lk< r 3 /■» _ r 4 

x^^' k (vKAi^' l (vWAl(r ] )-d^(r ] )]\U} + hM. . (14) 

P 
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As one can see, the dependence on the strong coupling constant disappears and the 
instanton contribution to the Sivers function is proportional to the packing fraction of 
instantons in the vacuum. We now use Eq. (5) for the quark fields; by taking also the 
Fourier transforms of the gluon fields, using the perturbative gluonic propagator (given 
consistently in a covariant gauge, see comments after Eq. (2) ) as an approximation for 
the correlator of two gauge fields, 

(A+mup)) = 5^( q+P ) (^^) , (is) 

and using Eq. (12) where appropriate, we get 



fi „. , MpEpit 2 f f ,!'», 1 /• d : '' k. /• <r-k. 



f 1T a (x,k ± )j = —— J — G (qp e , ^ ^ + (27r)3y (27r)3 

x 2n5{q°)5{k + - q + - k+)5 {2 \k ± - q ± - k 1± ) 

x EEE tIM tx I <i mi < m2 <l 3 <J T. ) (16) 

rrii p aijkl 

x 0L 1 (fci)7 O 7 + 0m 2 (fci + g)0L 3 ( fc 3)7 O (^7 + -7 + ^)0 m 4(fc3-g) + /i-c. . 

Concerning the hermitian conjugate contribution, one can easily see that it can be 
obtained by simply exchanging q <-> — q in the first term of the previous equation. Taking 
also into account the relation 

1 - = l (2n)5(q + ) , (17) 



q + — it q + + le 

performing the integrations over coordinates, and making use of the Dirac delta functions, 
one can simplify the sum of the two contributions, which finally reads 

f ±a, h \ .MpEp^f r f d"q ± G(q ±Pc ) f d 3 k, f d 3 k 3 1 

h T {x,k ± )j - t ^—— V2 J — _ 

x S (k + - kf - q + ) 5 i2 \k L - k 1± - q ± ) Tf m} 

rrii 

x 0j ni (fci)7V0m 2 (fc 1 + q ± )0^(fc 3 )r( g± )0 m4 (fc3-g ± ) , (18) 

where q± = \q±\, 

T { m} = EE W tx I <\ m y a>m2 a*l : /^ mi \ u > (19) 

/3 aijkl 

= 2 ^ ? tij^kli (S+z\ a% a,m 1 ^a,m 2 a l3]m 3 a p,rn i \S-z) + (S-z\ a a, mi a i,m 2 °/3,m 3 a /3,m 4 I ^+2) ) > 
(3 aijkl 

r(9x)=7°(7-^7 + -7 + 7-^), (20) 

and where a (not summed), (5 — u,d are flavour indexes and {m} stands for the set of 
quark helicity indexes (mi, m,2, m 3 , m 4 ). 

It is useful to compare this result with the perturbative contribution, Ref. [26]: 

f#(x, k ± ) pQCD = i 2g] ^ E Tf m} I^Af s p m (qD ^ 1>n «(<I-L, *0 > (21) 



where the current quark-gluon vertex contribution (common to both the perturbative and 
the instanton approach) is given by 

K u m 2 (Q±, *) = <PU k ~ q±h°l + <Pm 2 (k) , (22) 
and consists of two terms, with and without helicity flip of the struck quark: 

pc / i \ _ _J_ / rpc;even r , rpc;odd r 



v/2 

^A> 2 
\/2 



t (k)t (k') + ^t 1 (k)t (k>) + ^ to (k)t 1 (k') 



+ j^j [k ■ k' + i2m 1 (k x q ±y - k y q ±x )]t 1 (k)t 1 (k') 



+ 



+ ^^i,WMi')]t,-,l (23) 



c; odd 
mi 



where iV 2 = 4vriV 2 J R[>, fc = | fe|, k! = k-q ± , k' =\k'\. Notice that both F^ ven and 
have a symmetric, -F s c ym , and an antisymmetric, F£ ntisym , component for mi <-> —mi. 
The (perturbative) spectator quark-gluon vertex contribution reads 

/' d^k 1 
^ 0L 3 ( fc 3)7°7 + m4 (fc3 - q±) = ^= F s P {q\) 5 ms , m4 (24) 

where the spectator function FJ,(g 2 ) is given by (see also Ref. [26]) 

F s P {q\) = -^—^ J K 2 dKd(cos6) |t (K)*o(K) + ^— ^t (K)ti(/c') 

+ cosgt 1 (/t)t (/t / ) + - - q±cos6 t^t^K) ) , (25) 

and we have defined n = fc 3 , k = | fe 3 |, «/ = fc 3 — q ± , «/ = | k'| = [/t 2 +g 2 — 2Kg_ L cos#] 1 / 2 . 

It is worth to emphasize at this point that, according to our definition of the matrix 
elements (\ x \ . . . | | x ) (see comments after Eq. (4) and Eq. (19) ) the only non-vanishing 
T m 1 ,m 2 ,m 3 ,m 4 are those of the form T« i mm _ m3 or T^ _ mx mz mr Therefore, one-gluon 
exchange contributions, associated either with the perturbative or the nonperturbative 
spectator vertex, can be non-vanishing only when we have just (and at least) one helicity- 
flip term, either in the struck quark vertex, or in the spectator one, but not in both 
of them. In the perturbative case, the struck-quark helicity-flip situation takes place, 
Ref. [26] . As we will see below, the specific instanton- induced spin-spin correlation allows 
both the struck- and spectator- quark helicity flip mechanisms to contribute. 

In the nonperturbative case, it is convenient to define an instanton analogue of the 
function Fp(q\): 

F Ims , m4 (q±) = /7^^ 3 ( fc 3)r(g ± )0 m4 (fc 3 -g ± ) 

= 7y| ( Fl'%T(q±) <Wm 4 + Fl'ms(q±) <W-m 4 ) ' ( 26 ) 
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so that the instanton contribution to the Sivers function reads 
AV(^), = *^02 / ^±. G (q ±Pc ) T,Tl m} F^ m2 (q ± ,k) F/ m3jm4 ( 9± ). 

(27) 

Working in the standard representation of the Dirac matrices 



a 



-<= u o ) • ^ + - 7i i * *j ' 777 = 7i u; .3 j ' (28) 

one finds 

F^ m3 , m M±) =V2N 2 j ^Xl 3 {i[q± * <t] z t ( K )t ( K ') - (<r-<7 ± )(™') 1 t («)*i(«') 

+ (cr-K)(cr-qf ± ) - t 1 (/t)t (/t / ) - 2(<t-k)[<7_l x <r] z (<T-/« / ) *i(k)*i(k') } Xm 4 . 

(29) 

After some algebra and taking into account symmetry properties of the integrand, one 
gets 

F i'm 3 (<?_!_) = J2^y, J K dn d i cos v) Q± I cos 9 t {K )ti{K) — i («)*i(« ) I (30) 

2iV 2 /■ 
F i''m 3 d (q±) = (2m 3 g ±x - ig ±w ) y /t 2 ck d(cos#) 

x ^ (^o(«0 + K ~^ co&9 *i(k)*i(«')) • (31) 

Notice that Fj'™^ n is in fact independent of m 3 , while Fj'^ has both a symmetric 
and an antisymmetric component, r I sym and t 1 T antisym , * OT m 3 ^ — m 3- 

This time we can have two additive contributions to the Sivers function: Ff.°^ n and 
F s',odT ' i corresponding respectively to helicity-flip of the struck quark and no helicity flip 
for the spectator quark, and viceversa, 

fia/ K \ \ 1 / rpc;even ps; odd rpa , rpc;odd rps; even rpa \ /oo\ 

J IT V > "'-L/I / j ^ mi m 3 mi,mi,m 3 ,—m 3 ' mi m 3 mi,—mi,m 3 ,m 3 J ' \ ) 

mi ,m 3 

Therefore both the spectator and the current quark spin-flip terms contribute to the 
instanton- induced Sivers function. This is in contrast to the purely perturbative result, 
Eq.s (23) and (24), where only the current quark spin-flip gives a non- vanishing contri- 
bution to the Sivers function. Moreover, as we will see in a moment, it turns out that the 
instanton contribution has a completely different flavour dependence, due to the struc- 
ture of the matrix element in Eq. (27). This could open a new way of understanding the 
unexpected flavor dependence of the Sivers function, as extracted from recent HERMES 
data [2] (see below). 

An explicit calculation of the non-vanishing proton matrix elements T" m y in Eq. (19) 
shows that they have the same value for any allowed combination of mi, m 3 , so that, e.g., 

ppc; even jpc; even j-is; odd rpa _ rpa ST~^ pc; even rps; odd 

s;odd / j mi I m 3 mi,mi,m 3 ,—m 3 TTTI ' i mi I m 3 

mi,m 3 mi,m 3 

i i 

E T a ( _ \ "> ppc\even \ ( _ \ 1 rps; odd \ 

1 mi,mi,m 3 ,-m 3 \ r> mi / V 9 ^ m 'A ) 

mi,m 3 ^ mi ^ m 3 

/~ia\even rpc;even rps;odd fQQ^ 

— ^ilodd r sym r I symi 



where F^ en and Fj'" 1 ^ can be obtained from Eq.s (23) and (31) respectively. Analo- 
gously, 

rpc;odd /~ta\odd rpc;odd rps;even f^/H 

s;even I I even sym I ' V / 

where again F^ d can be obtained from Eq. (23). 

Notice that Cf\ °J^ n corresponds to the perturbative case, see Eq.s (21), (23), (24) and 
Ref. [26]. Straightforward calculation gives 



C7 



u I even 
odd 



"9 ' °" 



d\ even 
odd 



and 



u I odd 



16 , 



d\ odd 

even 



4 

+ 9- 



(35) 



(36) 



Therefore, in the instanton approach the total expression for the Sivers function has 
two contributions of a very different origin, which sum up together: 



x 



k ± )j = M P Ep ^N 4 j^y 5 J d 2 q ± ^ G(q ± p c ) 
to(k)t (k') + ^ tl (k)to(k') + ^t (k)t 1 (k') + k2 ~ ^ cos< 
n 2 dn d(cos9) 



even 
odd 



x | C? 

q± cos 



kk> 



t {K)t (K') + hiKjhin') 



K 



■h(k)ti(k') 
(37) 



+ C7 7 ° 



a | odd 

even 



k 



h(k)t (k') 



k L -q L cos(p , k z q± cos0 , 
— t {k)ti{k ) H — ti(fc)ti(fc ) 



fc' 



9± 



x — / k d,Kd(cos6) 
k± 



n , l\ , \ K COS & ~ Q± / x / /x 

COSUt (K )ti{K) to{K)ti{K ) 



where is the azimuthal angle of q ± , and 
uj(4x - 1) 



k z 



Ro 



(38) 



Notice that to go from Eq.s (23), (30), (31), (27) to Eq. (37) we have used the identity 



J d 2 q ± q ±i F(k 2 ,q 2 ± ,k ± -q ± ) = -± J d 2 q ± k ± q ± F(k 2 ,q 2 ± ,k ± -q ± ] 



(39) 



where F is any function with the functional dependence shown and % — x,y. In this way, 
we can effectively take in Eq.s (23), (30), (31), q±_i — > h L (k±-q ± )/k 2 L , and the numerator 
of the Sivers function, Eq. (27), becomes proportional to k y , which simplifies, as it must 
be, with the overall factor l/k y originally present in our expression. 



4 Numerical results 

In order to estimate the Sivers function for u and d valence quarks within the instanton 
approach proposed in this paper, we do not make any attempt to fix the parameters 
of the instanton model by fitting available results on single spin asymmetries. Rather, 
we consider widely accepted values for these parameters, coming from a comprehensive 
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Figure 2: The instanton (thick lines) and 
the one-gluon (thin lines) contributions to 
the first moment of the <i-quark (solid) and 
w-quark (dashed) Sivers function, vs. x. 
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k ± (GeV/c) 

Figure 3: The instanton and one-gluon 
contributions to the Sivers function at 
2=0.3 vs. fc^(GeV/c). Notations are the 
same as in Fig. 2. 



application of the instanton model to nonperturbative QCD (see Ref. [17] for more detail): 
p c ~ 0.3 fm, m* « 170 MeV and n c « 0.5 fm -4 . For the perturbative contribution we will 
take the following value of the strong coupling constant, a s ~ 0.3. This value is usually 
adopted for estimation of the one-gluon exchange contribution to SSA for HERMES 
kinematics (see e.g. Ref. [28]). This has to been taken into account when comparing 
our results with those of Yuan [26]. 

Using these parameters, in Fig. 2 we present the instanton and the perturbative, one- 
gluon exchange, contributions to the first moment of the Sivers function, 



/^(x) = J d'k ± k ± ) , (40) 

for u and d quarks, as a function of x. In Fig. 3, the Sivers function is shown as a function 
of k±, at x = 0.3. The main features of these results are that the instanton contribution 
is sizable and, in contrast with the purely perturbative case, it has the same (negative) 
sign and approximately the same magnitude for u and d quarks. 

As a consequence, when summing the two independent perturbative and instanton 
terms to get the total contribution to the Sivers function, for the <i-quark the two terms 
almost cancel each other, leading to a small, negative Sivers function. On the contrary, 
in the case of the w-quark Sivers function, the instanton contribution sums up with the 
perturbative one. This can be seen in Fig. 4, where we also compare the total contribution 
to the first moment of the Sivers function with recent parameterizations obtained by 
fitting the results of the HERMES Collaboration for the single spin asymmetry in semi- 
inclusive DIS pion production [2], see first of Ref.s [16]. From these results, one can 
observe that the present estimate of the instanton contribution to the valence quark 
Sivers function compares reasonably well, for u quarks, with parameterizations leading 
to a good description of the HERMES data in the valence-quark region, x > 0.3, where 
one can believe in the MIT-bag based calculations. On the other hand, for d quarks our 
results are smaller than one could expect from data fitting and differ from those of other 
model calculations and available parameterizations. 
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Figure 4: The total contribution to the first moment of d (solid line) and u (dashed line) 
quark Sivers distributions as a function of x. The l-a band for the parameterizations 
extracted by fitting the HERMES data (see first of Ref.s [16]) is also shown. 

5 Conclusions 

In this paper, the instanton contribution to the valence quark Sivers function has been 
estimated, adopting the MIT bag model for the nucleon wave function. Our results show 
that this contribution can be sizable, since it is not suppressed by powers of the strong 
coupling constant. The specific flavour dependence of the instanton contribution, as 
compared to that of the perturbative one, leads to a large, negative total contribution to 
the w-quark Sivers function. On the other hand, it suppresses the total contribution in 
the (i-quark case. However, several points should be kept in mind here: 

a) We have used widely accepted values for the basic parameters of the instanton model. 
These parameters are known with relatively large uncertainties which reflect into the 
results of our calculations. As an example, using the Diakonov and Petrov instanton liquid 
model estimate for the quark anomalous chromomagnetic moment (see e.g. Ref. [18]) 
would lead to an almost twice larger instanton contribution. 

b) The value adopted for the strong coupling constant plays a role in the interplay between 
the purely perturbative and the instanton terms, since the first contribution is proportional 
to a s while the second one depends on it in a more involved and indirect way. Changing 
the effective value of a s , which is not well known a priori, would therefore modify in 
particular the d quark case. 

c) Due to the effective instanton size, the instanton contribution could in principle lead to 
a different k± dependence of the Sivers function as compared to the purely perturbative 
case. In our results, however, the two contributions show a very similar behaviour in this 
respect. This is because in this first simplified approach the effective kj_ dependence of 
the Sivers function is basically determined by the MIT-bag model wave functions rather 
than by specific instanton properties. 

d) There could be in principle additional instanton contributions to the Sivers function, 
coming from more complex Feynman diagrams, involving e.g. the spectator diquark as 
a whole. Moreover, these contributions might be flavour-dependent, leading to different 
changes in the u and d quark Sivers function. 
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A detailed treatment of all these points is beyond the scope of this paper, which was 
mainly intended to show that instanton effects may be relevant for the Sivers distribution. 
A deepest study of instanton contributions to the Sivers function, including additional 
diagrams and based on more elaborated models for the nucleon wave function is in progress 
[32]. 
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